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(54) Precede de traitement d'une image incluant une etape de chainage, et appareil d'imagerie 
medicale incluant des moyens pour mettre en oeuvre ce precede 



(57) Procede de traitement cfimage incluant I'acqui- 
sition (10) d'une image d'intensite (J) formee d'une ma- 
trice de points (A) munie de valeurs d'intensite non bi- 
naires, representant des objets filiformes, et incluant 
des etapes de detection des objets filiformes 
comprenant : 

une determination automatique (20) d'un champ de 
vecteurs d'intensitd (0) formant une image 
d'ETATS constitues par des caracteristiques de 
grandeurs et d'angles fl|^||, 6] des vecteurs d'inten- 



site aux diff erents points (A) de I'image d'intensite, 
une operation de chainage (30) des points d'un ob- 
jet filiforme individuel au moyen d'un flttrage auto- 
matique [Q(i)] de I'image des ETATS qui effectue 
une selection d'un champ dense et continu de vec- 
teurs d'intensite correspondant a des points de I'ob- 
jet filiforme. 

Application : Appareils d'imagerie medicale 
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Description 

L'invention conceme un precede de traitement 
d'image incluant I'acquisition d'une image d'intensite 
formee cfune matrice de points munie de vaieurs (fin- s 
tensite non binaires, representant des objets filiformes, 
et incluant des etapes de detection des objets filiformes. 

L'invention trouve une application dans I'industrie 
de fabrication des appareils de radbtogie aux rayons X 
par exemple. 

Une methode de detection des bonds est deja con- 
nuede la publication : "IEEETRANSACTIONSON PAT- 
TERN ANALYSIS AND MACHINE INTELLIGENCE, vol. 
14, N°9, SEPT 1992, Object Delineation in Noisy Ima- 
ges by a Modified Policy-Iteration Method, by Adrie C. 
M. Dumay, Harielle N.A.J. Claessens and alii, pp. 
952-958'. 

Cette publication decrit un algorithme de detection 
de ligne de contour d'objets isoles dans des images 
bruitees qui est une variante cfune methode bien con- 
nue de Phomme du metier sous le nom de programma- 
tion dynamique. La methode de programmation dyna- 
mique s'applique sous deux conditions : e'est une me- 
thode de recherche du chemin le plus court pour relier 
un groupe de points avec un cout minimal. Cette me- 
thode a besoin d'etre guidee ; elle n'est done pas auto- 
matisable. 

La presente invention a pour but de proposer un 
procede de traitement d'image, tel que defini plus haut, 
s'appliquant sous un plus grand nombre de conditions 
que la methode de programmation dynamique connue : 
les etapes de calcut doivent se derouler de maniere 
automat ique ; le nombre des etapes doit etre aussi fai- 
ble que possible ; ces etapes doivent determiner le che- 
min le plus court pour relier un groupe de points appar- 
tenant a un objet filiforme individuel, avec le coOt le plus 
faible possible ; et le chemin choisi doit inclure le plus 
grand nombre possible de points de I'objet filiforme, 
done etre le plus long et le plus regulier possible dans 
I'image. 

La condition selon laquelle le chemin trouve doit in- 
clure la plus grande densite possible de points apparte- 
nant a I'objet filiforme est inconnue de la methode de 
programmation dynamique. 

Selon l'invention, ce probleme est resolu par un pro- 
ceed de traitement d'image incluant Pacquisition cfune 
image d'intensite formee d'une matrice de points munie 
de vaieurs d'intensite non binaires, representant des ob- 
jets filiformes, et incluant des etapes de detection des 
objets filiformes comprenant : 

une determination automatique d'un champ de vec- 
teurs d'intensite (7) formant une image d'ETATS 
constitues par des caracteristiques de grandeurs et 
d'angles fl|^||,9] des vecteurs d'intensite aux drffe- 
rents points de I* image cfintensite, 
une operation de chatnage des points d'un objet fi- 
liforme individuel au moyen d'un filtrage automat i- 
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que [Q(i)J de i'image des ETATS qui effectue une 
selection cf un champ dense et continu de vecteurs 
d'intensite correspondant a des points de I'objet fi- 
liforme. 

Ce procede possede I'avantage de s'appliquer au 
traitement d'images medical es de rayons X, qui sont 
peu contrastees a Porigine, telles que des arteriogram- 
mes. Dans un arteriogramme pulmonaire, ce procede 
permet par exemple la segmentation automatique d'ar- 
teres de differents diametres par rapport au fond qui 
comprend une representation genante des cdtes. Dans 
un arteriogramme c6 rubral, ce procede permet par 
exemple une segmentation automatique cfarteres tres 
enchevetrees. D'une maniere generate ce procede off re 
I'avantage de permettre la detection et la reconstruction 
amelbrees des bords ou des lignes dans de nombreux 
types d'objets representee dans des images, sans con- 
dition d'origine de ces images. 

Un appareil d'imagerie medicale incluant des 
moyens pour mettre en oeuvre ce procede comprend : 

un systeme d'acquisrtion de donnees de represen- 
tation d'une image sous forme d'une matrice bidi- 
mensionnelle de points ayant des vaieurs d'intensi- 
tes non binaires, reperes par leurs coordonnees 
dans la matrice, 

un systeme de traitement d'image comprenant un 
microprocesseur ayant acces aux donnees de re- 
presentation de I'image et munis de moyens pour 
mettre en oeuvre les etapes du procede, 
et un systeme cfaffichage et/ou d'enregistrement 
des donnees d'image acquises et/ou traitees par le 
procede. 

Un proc6d6, et un appareil pour mettre en oeuvre 
ce procede sont decrits ci-apres en detail, i (lustres par 
les figures schetnatiques annexees parmi lesquelles : 

les FIGS 1 A et 1 B illustrent le procede par des dia- 

grammes de blocs fonctionnels, 

les FIGS 2A et 2B illustrent des d6terminations 

d'ETATS, 

la FIG. 3 illustre la discretisation de la caracteristi- 
que d'angle, dans la determination d'ETATS, 
les FIGS.4Aa 4D montrent un exemple de balayage 
automatique de I'image des ETATS, 
les FIGS.5A, 5B representent des secteurs de re- 
cherche, 

la FIG. 6 illustre un pistage par le procede, 

la FIG. 7 A montre un arteriogramme original, et la 

FIG.7B montre cet arteriogramme traite selon le 

procede pour chaTner les points de bords de vais- 

seaux, 

la FIG.8 montre un appareil a rayon X mettant en 
oeuvre le procede. 
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PROCEDE 

Le procede d6crit ci-aprfcs est il lustre par les FIGS. 
1 A et 1 B sous forme de diagrammes de blocs fonction- 
nels, et conceme d'une maniere generate le traitement 
d'une image d'intensite non binaire, appelee image ori- 
ginate J, telle que montree a tit re cf exemple sur la FIG. 
7A, representant des objets filiformes. Les objets filifor- 
mes peuvent etre, a titre d'exemple, des bords d'objets, 
des aretes, des lignes de centre cfobjets. Ce procede 
comprend des Stapes de traitement appliqu6es a cette 
image originate J et incluant au moins une etape de 
chatnage 30 de pixels appartenant a un objet filiform e. 
Le precede comprend en outre favorablement une 6ta- 
pe de pistage 40 de chaque objet filiforme tndividuelle- 
ment dans I'image K tel qu'iltustre par la FIG.7B. 

En reference a la FIG. 1 A, le precede comprend : 

1) une premiere etape preliminaire 10 d'acquisition 
de Timage originale J atraiter. Cette image originate 
J peut etre acquise sous forme d'une mat rice bidi- 
menstonnefle de pixels ou de points. En reference 
aux FIGS.3A, 3B, chaque pixel A est muni de coor- 
donnees x,y de localisation dans la matrice de Tima- 
ge J par rapport a des axes de coordonnees O)? , 
O? et est muni d'une valeur d'intensite non binaire 
l(x,y), qui peut etre mesuree ou calculee. Le prece- 
de n'est pas limite au traitement cfimages bidimen- 
sionnelles, mais peut aussi etre applique par exem- 
ple au traitement cfimages tridimensionnelles, ob- 
tenues par le moyen de systemes cfimagerie parti- 
culiers qui permettent de former des images de vo- 
lumes. L'homme du metier est en mesure de gene- 
ral iser facilement le procede deer it a titre d'exemple 
pour une image bidimensionnelle a une image tri- 
dimensionnelle. 

2) une seconde etape preliminaire 20 de determi- 
nation en chaque pixel A de caracteristiques des- 
criptives de la valeur d'intensite de ce pixel A par 
rapport aux valeurs d'intensite des pixels de son en- 
vironnement. Ces caracteristiques appetees ci- 
apres "ETAT" comprennent les coordonnees de lo- 
calisation x,y du pixel A associees a des caracte- 
ristiques locales orientees d'intensite detenminees 
de maniere automatique. 

A cet effet en reference a la FIG. 2 A, par exemple 
pour la detection d'un objet filiforme qui est une ligne de 
bord cfun objet quelconque dans I'image originale J, les 
caracteristiques locales orientees d'intensite sont favo- 
rablement le module du gradient d'intensite [| 3D en A, 
et ('orientation <p du vecteur gradient (J calcuie par rap- 
port a Tun des axes de coordonnee . Ces calculs 
peuvent dtre eff ectues favorablement a titre d'exemple, 
par balayage conventionnel systematique de Timage 
originate, au moyen de filtres de derivation des vaieurs 
d'intensite appeles filtres de Sovel qui determinent la 
pente notee Gx des vaieurs d'intensite a partir cfun pixel 



donn6 A dans la direction Ox* , et la pente notee Gy des 
vaieurs d'intensite a partir du meme pixel A dans la di- 
rection U? . 

Les composantes du vecteur gradient, ou pontes 
5 d'intensite, Gx et Gy foumissent, a Tissue de I'etape 20 : 
le . module du gradient d'intensite 

•l|3ll = 7 G x + Gy 

Tangle <p que fait le vecteur gradient (? avec Taxe 
O)? pris comme repere : 9 = Arctg (Gy/Gx) 

10 On dispose ainsi d M ETATS" relatifs aux pixels de 
Timage d'origine J, qui torment un "champ de caracte- 
ristiques locales orientees d'intensite", ou champ 
d'ETATS, ou champ de vecteurs d'intensite formant une 
image d'ETATS. 

is Dans le cas ou ('acquisition d' image a ete faite en 
trois dimensions, alors le gradient est calcuie dans un 
volume au lieu d'etre calcuie sur un plan, et la suite du 
traitement est adapte a ces donn6es de toute maniere 
automatique connue de Thomme du metier. 

20 L'invention n'est pas Iimit6e a la detection des li- 
gnes de bord. Toute ligne peut etre d6tect6e de maniere 
automatique pourvu qu'on associe a chaque localisation 
de pixels dans Timage un vecteur d'intensite, detini par 
son module et son angle d'orientatbn, et qui est des- 

25 criptif de la ligne dans son environnement. 

En reference a la FIG.2B, dans un second exemple, 
on peut chercher a detecter des lignes metJianes de 
vaisseaux sanguins dans un arteriogramme. A cet effet, 
on peut par exemple d'abord determiner les lignes de 

30 bords des vaisseaux par une methode connue des gra- 
dients. On peut ainsi determiner deux points E1 , E2 de 
bords opposes B1 , B2 d'un vaisseaux dont on a calcuie 
respectivement le vecteur gradient (a E1 associe a un an- 
gle q> E1 , et le vecteur gradient C?^ associe a un angle 

35 9^. Un point de la ligne mediane du vaisseaux est alors 
A(x,y) situe au point d'intersection des directions des 
vecteurs gradient fj E1 , (S^. Ce point est muni de ca- 
racteristiques locales orientees constituees par son in- 
tensite l(x,y) calculee ou mesuree, et par un angle y que 

40 fait la bissectrice des vecteurs gradients (? E1 et (? £2 
avec le rep6re O)? . On d6finit ainsi dans Tetape 20 pour 
les points medians potentiels, une image d'ETATS for- 
m6e d'un "champ de caracteristiques locales orientees 
d'intensite" ou champ d'ETATS ou champ de vecteurs 

45 d'intensite, 

En reference a la FIG. 3, on peut choisir pour sim- 
plifier les calculs de discretiser Tangle 9 ou \\t du vecteur 
d'intensite consid6r6, par exemple selon des directions 
0, a 0 Q et + n a e 8 + tc, de n/8 en tc/8. 

50 D'une maniere generate, un ETAT est defini au 
moins par une localisation x,y et un vecteur d'intensite 
7 dont on utilise deux composantes qui sont le module 
d'intensite du vecteur et un angle G cforientation de 
$ par rapport a un repere Ox* , de preference discretise. 

55 3) Une etape 30 de chatnage des pixels apparte- 
nant a un objet filiforme de Timage d'origine. Le chatna- 
ge est une operation qui consiste a determiner a partir 
d'un pixel donne, de proche en proche, un chemin forme 
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de pixels rfun objet filifonme qui est le plus long et le plus 
regulier possible dans Timage d'origine. A cet effet, on 
determine cfabord, pour un pixel donne, un pixel pred6- 
cesseur (ou antecedent) qui appartient a un tel chemin 
le plus long et le plus regulier possible, puis on determi- 
ne a partir de ce predecesseur un predecesseur du pre- 
decesseur appartenant a ce chemin. On associe ensui- 
te de proche en proche tous les pixels ayant des ETATS 
qui peuvent dtre trouv6s avec cette condition, et on ob- 
tient une chain e de pixels representant au mieux I'objet 
filiforme correspondant dans Timage d'origine. On don- 
ne ci-apres un critere de selection qui permet de deter- 
miner, a partir d'un ETAT lie a un pixel donne, I'ETAT 
predecesseur, s'il en existe, lie a un pixel predecesseur 
qui fait partie d'un tel chemin le plus long et te plus re- 
gulier possible. A partir cfun ETAT lie a un pixel donne, 
ce critere 6 limine les ETATS Ii6s a des pixels pred6ces- 
seurs potentiels qui conduiraient a la construction de 
chemins courts et/ou irrdguliers, et determine, s'il en 
existe, un ETAT lie a un pixel predecesseur potentiel qui 
permet de construire le chemin le plus long et regulier 
possible. 

Selon I'invention, vis-a-vis de T6tat de la technique, 
dans le but d'ameiiorer la detection des objets filiformes, 
on a ainsi pos6 des contraintes plus severes d'apparie- 
ment d'un ETAT donne et d'un ETAT antecedent, tout 
en imposant que ('operation de chaTnage 30 se d6roule 
au moyen d'un nombre de sous-stapes aussi faible que 
possible, de maniere automatique. Les problemes sont 
resolus par la mise en oeuvre d'un fittrage de Timage 
d'ETATS avec un filtre applique en un pixel donn6 ayant 
un ETAT i, ce filtre 6tant deTini par un noyau exprim6 
selon la formule survante : 

xpPGJexpC-CjjJ.ouC^La^ (1a) 

d'ou il resulte un critere de selection d'un ETAT prede- 
cesseur potentiel d'un ETAT donn6 exprim6 par : 

QC^ZpPOJexp^) (1b) 

ou Q(i) est la valeur cf intensite filtree cfun ETAT t relatif 
au pixel donne, par exemple la valeur de || v*|| en i apres 
fittrage not£e ||7j| F ; ou P(j) est la valeur d'intensite non 
filtree d'un ETAT j relatif a un pixel antecedent potentiel, 
par exemple la valeur non filtr6e de ||^|| en j, notee ; 
ou a, m est une mesure de disparite entre des ETATS 
adjacents, note l,m, situes sur le chemin reliant i a j ; ou 
Cjj est le cout du chemin reliant i a j obtenu en evaluant 
une somme Z des disparites entre les ETATS adjacents 
situ6s sur le chemin reliant i a j. Les chemins retenus 
sont determines par Toptimisation du critere (1) qui est 
obtenu e lorsque le cout du chemin est petit et la den- 
site des ETATS sur ce chemin est grande. L'optimisation 
de ce critere permet de determiner ledit chemin le plus 



long et le plus regulier. 

II apparaTt que la mise en oeuvre du critere (1) est 
trop consommatrice en temps de calcul. C'est pourquoi 
le precede selon Tinventkxi propose la mise en oeuvre 
5 du critere recursif suivant : 

Q(i)=pP(i) + aQ(j) (2) 

10 qui sur le plan mathematique est completement equiva- 
lent au critere (1 ), comme peut le deduire facilement 
Thomme du metier. Le fi It rage au moyen du critere (2) 
est un algorithme recursif qui effectue une somme de P 
(i) qui est TETAT i non filtre par exemple ||tf J|, et de Q(j) 

75 qui est un ETAT filtre lie a un pixel j predecesseur retenu 
note || v*j|| F pondere par un facteur de probabilite cfappa- 
riement a = cty pour foumir Q(i) qui est TETAT i filtre par 
exemple |(v* J| F . Le critere 2) est mis en oeuvre au moyen 
d'un nombre d'operations reduit vis-a-vis du critere 1) 

20 comprenant, en reference a la FIG . 1 B des sous-etapes 
de: 

definition 31, pour deux ETATS donnes i et j, d'une 
probabilite d'appariement ct^ que Ton 6crit : ot^ = exp 
25 (- a { j) ou a^ est une mesure de disparite entre deux 
ETATS i,j donnes ; 

initialisation 32 du critere Q(i) hxalement selon la 
formule Q(i) = p.P(i) ou Q(i) est la valeur d'intensite 
de TETAT i filtre, not6e ||$ J| F , ou p est une constante 
30 et ou P(i) est la valeur d'intensite de TETAT i non 
encore filtre, notee ||$J|, . 
balayage 33 de Timage des ETATS, d'une maniere 
conventionnelle selon au moins une direction notee 

zt. 

35 recherche 34, dans la zone d6ja balay6e, d'un ETAT 
predecesseur, note ETAT k(i), pour lequel le critere 
Q(i) est maximal, 

stockage de TETAT k(i) correspondant a ce maxi- 
mum. 

40 

Le filtrage selon le critere Q(i) detecte des ETATS 
continus, e'est-a-dire un champ de vecteurs continu. 
Par exemple dans le cas de la detection des bords, ce 
filtrage detecte les champs de gradients continus. Ce 

45 filtrage determine TETAT predecesseur k(i) le plus pro- 
bable qui doit §tre lie a TETAT i, de maniere automatique 
et a moindre cout. 

En reference a la FIG. 4 A, pour rendre retape de 
chaTnage automatique, Timage des ETATS est baiayee 

so automatiquement seton au moins une direction. Par 
exemple un balayage conventionnel peut etre effectue 
selon les lignes dans la direction ZT paralieiement a 
Taxe O? , du haut a gauche vers le bas a droite, dans 
le sens CAUSAL, dans Timage des ETATS. Ce balayage 

55 permet de determiner TETAT predecesseur k(i) dans la 
region de Timage deja baiayee, qui se situe plus haut 
ou plus a gauche que TETAT courant i bataye. 

De preference en reference a la FIG.4B, apres un 
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premier passage CAUSAL selon Zt parailele a O? , 
i'image des ETATS est ensuite balay6e automat ique- 
ment selon les lignes dans la direction Z? ANTICAUSA- 
LE, du bas a drorte vers le haul a gauche de I'image des 
ETATS, parallelement a O? . Ce second balayage per- 
met de determiner un ETAT h(i) successeur de PETAT 
courant i balaye. 

De preference, en reference aux FIGS 4C et 4D, 
des balayages systematiques causa ux et anticausaux 
sont encore effect u 6s selon des directions Z3 et 7Z per- 
pendiculaires aux directions ZT et Z? des premiers ba- 
layages. Dans I'exemple cite ou les premiers balayages 
systematiques sont farts selon les lignes paralleles a 
I'axe OX* , on effectue les seconds balayages systema- 
tiques selon des colonnes paralieies a I'axe oT , de ma- 
ni6re CAUSALE ET ANT1CAUSALE. 

En reference a la FIG. 2, une des caracteristiques 
du champ d'ETAT etant un angle 9 d'un vecteur v* re- 
p6r6 par rapport a I'axe 0)t par exemple, on determine 
un voisinage dans la region d'image d'ETATS d6ja ba- 
layes qui est constitue par un secteur angulaire de re- 
cherche note NB centre sur retat i balaye et qui fait un 
angle ±A0autourde la direction 0du vecteur^ . Grace 
a ce voisinage NB, on peut rechercher le meilleur ETAT 
voisin dans un champ de vecteurs continu. Dans le cas 
de detection des bords, on cherche un ETAT j qui maxi- 
mise le critere 2), c'est-a-dire dont I'angle 6 est proche 
de celui de I'ETAT i donn6, et pour lequel la mesure de 
disparite a^en cooperation avec la valeur du module du 
gradient ||Gj|p en j et la valeur du module du gradient 

en i permettent d'obtenir || 0? J| F maximal en appli- 
quant le critere 2). 

On peut favorablement choisir un secteur angulaire 
deiimite par ±A0 = ±rc/8, lorsque Pangle 0 est discr6tis6 
de it/8 en Selon ('orientation 0 du vecteur v*j a 
I'ETAT courant i, on dispose cfautant de secteurs de re- 
cherche NB que d'orientations 0 discr6tis6es telles que 
montr6es sur la FIG. 3. 

Dans le secteur de recherche NB, pour que la me- 
sure de disparite a^ reste raisonnable, le voisinage est 
egalement limite par une distance maximale Rde I'ETAT 
voisin examine j vis-a-vis de I'ETAT courant i balaye. II 
peut, par exemple, dans ce cas, se trouver 10 a 100 
ETATS j voisins possible dans le secteur angulaire NB 
determine. A cet effet, pour chaque ETAT courant i ba- 
laye, 10 a 100 essais correspondants sont effectues, 
dans la region d6ja balay6e et sur les ETATS j contenus 
dans le voisinage determine NB en appliquant le critere 
2). Parmi tous les ETATS j essayes dans le voisinage 
de recherche NB on retient celui qui maximise le critere 
Q(i), et on le note k(i). 

D'une maniere g£n6rale le chaTnage consiste a 
trouver, pour un ETAT courant i, dans une region de 
I'image des ETATS d6ja balayee, un ETAT k(i) parmi les 
ETATS j d'un secteur de recherche NB dont la direction 
0 est determinee en dependance de I'angle 0 lie a 
I'ETAT i courant, tel que les caracteristiques locales 
orientees liees a k(i) permettent d'optimiser le critere 2) 



selon ('expression : 

ll% = Pll^l + a ij ||^ j D F 3) 

5 

II en r6sulte que le meilleur ETAT voisin k(i) de 
I'ETAT donne i, dans ce voisinage, selon ce balayage et 
par ^application de ce critere, est localise a une distance 
interieure a R et est caract6ris6 par un vecteur d'inten- 

10 site dans un champ continu de vecteurs et est appeie 
pr6d6cesseur k(i) pour dtre chain6 a I'ETAT i. On stocke, 
dans la sous-6tape 35, 1'information que cet ETAT k(i), 
qui produit Q(i) fiftr6 maximal pour I'ETAT i, est le 
meilleur antecedent ou predecesseur de I'ETAT i. 

15 Apres les sous-6tapes 33 a 35 commencant par un 
balayage CAUSAL 33, on effectue une repetition des 
sous-etapes 33 a 35 en commencant par un balayage 
ANTI CAUSAL, ce qui permet de determiner par la m6- 
me methode un ETAT successeur qui est stocke et note 

20 h(i). On effectue ensuite des repetitions des sous-eta- 
pes 33 a 35 avec des balayages causaux et anticausaux 
perpendiculaires, et on cherche s'il existe un ETAT j, tou- 
jours dans le voisinage NB predetermine, qui foumirait 
une valeur Q(i) ameiioree. Dans ce cas, c'est cette nou- 

25 velle valeur cf ETAT j qui est stock6e pour k(i) ou h(i), et 
la valeur trouv6e dans les sous-etapes 33 a 35 avec des 
balayages precedants est oubliee. 

L'invention n'est pas limitee a apparier des ETATS 
determines a la m§me resolution. On peut favorable- 

30 ment filtrer i'image originate a diverses resolutions, ap- 
peiees echelles a, pour former un ensemble d'images 
multiechelles, et dans ce cas, les meilleurs ETATS voi- 
sins sont recherches aux differentes echelles, les echel- 
les les plus proches 6tant neanmoins priviiegi6es, par 

35 exemple en appliquant sur la valeur Q(j) dans l'6quation 
2), un poids approprie. 

Le chemin determine par la methode de chainage 
proposee possede des proprietes specifiques ; c'est 
pourquoi, on a donn6 au filtrage qui permet d'obtenir ce 

40 chemin le nom de FILTRAGE GEODESIQUE par ana- 
logie a une notation geographique. La determination de 
ce chemin permet de combler des petrtes coupures de 
I'objet filiforme dans I'image du fait que le meilleur ETAT 
voisin trouve dans le voisinage determine n'est pas obli- 

4S gatoirement I'ETAT qui est a la distance la plus courte, 
mais I'ETAT qui est a une distance courte en dependan- 
ce de caracteristiques locales orientees d'intensite de 
faible disparite. En revanche, le precede ne permet pas 
de combler des grandes coupures ce qui permet de s6- 

50 parer deux objets filiformes diff6rents. Ce filtrage GEO- 
DESIQUE est •contrainf par les conditions de localisa- 
tion, d'angle, et d'intensite et occasionne lie ment 
d'echelle qui sont liees au balayage, au secteur de re- 
cherche, et au critere 2) propose, et rehausse les objets 

55 allonges ou filiformes. 

4) Une etape 40 de pistage (back-tracking) automa- 
tique. Dans cette etape, on chaTne les ETATS corres- 
pondant aux differents objets filiformes de I'image. En 
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reference a la FIQ.6 le pistage automatique comprend 
une premiere determination d'un point i de I'image pour 
lequel on a trouve la vateur Q(i) la plus grande. On cree 
d'une part une chaine amont AM avec le predecesseur, 
puis le predecesseur du predecesseur etc... et on opere 
de proche en proche jusqu'a ce que i'on ne trouve plus 
de predecesseur. Et d'autre part, a partir de ce meme 
point i, on cree une chaTne aval AV avec le successeur, 
le successeur du successeur etc... et on opere de pro- 
che en proche jusqu'a ce qu'on ne trouve plus de suc- 
cesseur. Puis le pistage automatique comprend une se- 
conde determination d'un autre point de I'image pour le- 
quel on a trouve une valeur Q(i) aussi grande que pos- 
sible, a partir duquel on cree une nouvelle chaine amont 
et une nouvelle chaTne aval. Des pistages automat iques 
sont farts tant que Ton trouve des points de depart de 
chatnage ayantdes valeursQ(i) suffisamment grandes : 
on peut se fixer un seuil. Des chemins d6ja visites ne 
sont pas revisites. A cet effet, les nouveaux points de 
depart de chatnage possible ne sont retenus que s'ils 
sont a une distance suffisante pred6termin6e des chaT- 
nes d6ja existantes. Ainsi, des objets filiformes qui se 
croisent sont detectes indrvidueilement dans I'image. II 
ne peut y avoir d'ambiguite dans la segmentation des 
objets. Dans une 6tape de visualisation 50, la FIG.4B 
montre dans une image K le resurtat du chatnage des 
objets filiformes de I'image J de la FIG.4A. 

DISPOSITIF 

En reference a la FIG . 8, a tit re d'exemple, un appa- 
reil de radbgraphie numerisee comprend une source 1 
de rayons X; une table mobile 2 pour recevoir un 
patient ; un dispositif intensificateur d'image 3 couple a 
un tube video 4, qui foumit des donn6es a un systfcme 
8 de traitement num6rique d'image comprenant un mi- 
croprocesseur 5. Ce dernier comprend plusieurs sorties 
dont une sortie 6 coupiee a un monrteur 7 pour la visua- 
lisation de i'image radiographique ou des images d'in- 
tensite. 

L'image radiographique numerisee peut contenir 
512x512ou 1024x1024 pixels codes sur8bits ou 10 
bits. Chaque pixel peut etre ainsi affecte d'un parmi 256 
ou 1024 niveaux d'intensite. Par exemple, les regions 
sombres ont un niveau d'intensite faible, et les regions 
claires de I'image ont un niveau d'intensite 6lev6. L'ima- 
ge numerisee peut etre obtenue en mode fluoroscopi- 
que. (.'invention peut trouver son application particulie- 
rement dans le traitement d'images art6riographiques. 

Le traitement de I'image numerique selon le prece- 
de decrit plus haut est mis en oeuvre dans ses diff eren- 
tes etapes et operations dans le systeme 5. Les don- 
n6es sont stock6es dans une zone m6moire non repre- 
sentee et afflch6es si besoin sur un moniteur 7. Des 
moyens d'enregist rement 6galement notes 7 peuvent 
etre en outre utilises. La num6risation de rimage a trai- 
tor apporte un avantage dans la mise en oeuvre du pre- 
sent precede. 



Revendications 

1 . Precede de traitement d'image incluant I'acquisition 
(10) d'une image d'intensite (J) form6e d'une matri- 
5 ce de points (A) munie de valeurs d'intensite non 
binaires, representant des objets filiformes, et in- 
cluant des etapes de detection des objets filiformes 
comprenant : 
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15 



20 



25 



30 



35 



40 



une determination automatique (20) d'un 
champ de vecteurs d'intensite (v*) formant une 
image d'ETATS constitu6s par des caracteris- 
tiques de grandeurs et d'angles DI^O.G] des vec- 
teurs d'intensite aux differents points (A) de 
i'image d'intensite, 

une operation de chatnage (30) des points d'un 
objet filiforme individuel au moyen d'un filtrage 
automatique [Q(i)J de I'image des ETATS qui ef- 
fectue une selection d'un champ dense et con- 
tinu de vecteurs d'intensite correspondant a 
des points de i'objet filiforme. 

2. Procede selon la revendication 1 selon lequel Cope- 
ration de chatnage (30) par filtrage automatique de 
I'image des ETATS comprend : 

la definition d'un critere [Q(i)J de filtrage qui 
foumit la valeur fittreo de la grandeur []|v*J| F ] du 
vecteur d'intensite en un point d'ETAT(i) donn6 
de I'image des ETATS comme une somme de 
la grandeur [P(i),|| v* J|] non filtr6e du vecteur d'in- 
tensite en ce point et de la grandeur [QQ),||v^| F ] 
filtr6e du vecteur d'intensite en un point d'ETAT 
(j) voisin pond6ree par un coefficient de proba- 
bilite (ccy) d'appariement des points d'ETAT (i, 
j) en dependance d'une mesure de leur dispa- 
rite (a^), 

la recherche d'un point [k(i) ; h(i)] parmi des 
points d'ETAT (j) voisins dont les caract6risti- 
ques de grandeur et d'angle des vecteurs d'in- 
tensite maximise le critere [Q(i)] de filtrage au 
point d'ETAT donn6 (i). 



3. Precede selon la revendication 2, selon lequel 
45 ('operation de chatnage comprend : un balayage 

(33) systematique et automatique de I'image des 
ETATS selon au moins une direction (ZT) ; en un 
point d'ETAT (i) courant, une recherche (34) dans 
la zone de I'image d6ja balayee, d'un point d'ETAT 
so predecesseur [k(i)] et/ou successeur [h(i)] pour le- 
quel le critere [Q(i)] est maximal ; et le stockage 
desdits points d'ETAT predecesseur [k(i)j et/ou suc- 
cesseur [h(i)J comme points a chatner avec ledrt 
point courant (i). 

55 

4. Procede selon la revendication 3, selon lequel, 
('operation de recherche (34) d'un ETAT predeces- 
seur [k(i)] et/ou successeur [h(i)] d'un ETAT donn6 
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(i) dans I'image des ETATS comprend la determina- 
tion d'un voisinage de recherche (NB) forme d'un 
secteur angulaire centr6 sur ledit point d'ETAT don- 
ne (i), et qui couvre un angle donn6 (± AG) autour 
de Tangle d'orientation (0) du vecteur d'intensit6 5 

relatif au point d'ETAT donne (i), et la ddtermi- 
nation dans ce voisinage de recherche, d'un point 
d'ETAT [k{i) ; h(i)] parmi les points d'ETATSfj) du 
secteur de recherche pour lequel le critere est maxi- 
mal. io 

5. Procede selon la revendicatton 4, selon lequel le 
voisinage de recherche (NB) est limits en outre a 
une zone de rayon predetermine (R) autour du point 
courant (i). is 



6. Proced6 selon i'une des revendications 1 a 5, qui 
comprend en outre une operation de pistage (40) 
des points d'ETATS de drfferents objets filiformes 
indrviduellement. 20 



7. Procede selon la revendicatbn 6, selon lequel, 
Topgratbn de pistage comprend la determination 
d'un point d'ETAT (i) d'interfit, par Pexamen de la 
valeur [Q(i)J obtenue par la mise en oeuvre du cri- 2s 
tere de tlttrage, comme point de depart cf une ope- 
ration de chaTnage ; le chatnage de son precteces- 
seur [k(i)] et de son successeur [h(i)] et de proche 

en proche de tous les pr£decesseurs des pr&leces- 
seurs et successeurs des successeurs s'il s'en trou- 30 
ve, de part et d'autre (AM.AV) du point de depart, 
de maniere a former une chaine de points repr6- 
sentatifs de Tobjet filiforme ; le renouvellement de 
ces operations pour pister les d iff e rents objets fili- 
formes del'imageavec la condition que les point de 35 
depart sont choisis non confondus avec des points 
d'une chaine deja existante. 

8. Appareil d'imagerie m6dicale incluant des moyens 
pour mettre en oeuvre un proced6 selon I'une des 40 
revendications 1 a 7, comprenant un systeme d'ac- 
quisrtion (4) de donnees de representation d'une 
image sous forme cfune matrice bidimensionnelle 

de points ayant des valeurs d'intensites non binai- 
res, rep6res par leurs coordonnees dans la matrice, 45 



un systeme de traitement d'image (8) compre- 
nant un microprocesseur (5) ayant acces aux 
donn6es de repr6sentation de I'image et munis 
de moyens pour mettre en oeuvre les 6tapes so 
du precede, 

et un systeme d'afflchage et/ou d'enregistre- 
ment (7) des donnees d'image acquises et/ou 
traitees par le precede. 
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FIG. 1A 




FIG. 1B 
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FIG. 2B 
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